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ABSTRACT 

Current helicity quantifies the location of twisted and sheared non-potential structures in a magnetic field. 
We simulate the evolution of magnetic fields in the solar atmosphere in response to flux emergence and shearing 
by photospheric motions. In our global-scale simulation over many solar rotations the latitudinal distribution of 
current helicity develops a clear statistical pattern, matching the observed hemispheric sign at active latitudes. 
In agreement with observations there is significant scatter and intermixing of both signs of helicity, where 
we find local values of current helicity density that are much higher than those predicted by linear force-free 
extrapolations. Forthcoming full-disk vector magnetograms from Solar Dynamics Observatory will provide 
an ideal opportunity to test our theoretical results on the evolution and distribution of current helicity, both 
globally and in single active regions. 
Subject headings: Sun: activity — Sun: magnetic fields 



1. INTRODUCTION 

The twist of the solar magnetic field plays an important role 
in transient phenomena such as solar flares and coronal mass 
ejections, and in the dynamo processes that cause the 1 1 -year 
solar activity cycle. The magnetic twist can be measured in 
various ways. Magnetic helicity is an integral that quantifies 
topological complexity of field lines, such as linking, twist, or 
kinking (Bergejl99&[[999). For a closed magnetic system 
it is defined by H m = J A ■ Rd^x, and alternative definitions 
have been developed fo r open systems (Ber ger & Field[l984c 
lFinn&Antonsenlll985l) . 

In this letter we consider current helicity, which we define 



as 



B 2 ' 



(1) 



where B is the magnetic field and j = V x B is the current 
density. The quantity a has the advantage that it describes the 
local distribution of twist and shear in the magnetic field, and 
that it is more readily determined from limited observational 
data than H m which requires global information. For a force- 
free field (j x B = 0) we have j = aR and a, which may be a 
function of space, is a fundamental parameter that describes 
the torsion of the field lines around one another. Note that 
we shall not consider the integral current helicity H c = J j • 
Bd 3 x because unlike Hm it is not a near-conserved quantity in 
MHD (Demoulin 2007), and it does not even in general take 
the same sign as H m (except for linear force-free fields where 
a is constant in space and a, H c , and H m all have the same 
sign, Hagyard & Pevtsovll 19991) . 

There are two main techniques for estimating a from ob- 
served vector magnetograms, which so far only cover a small 
region of the solar surface such as a single active region: 



1. Compute j z = dBy/dx 
which should give a 



dB x /dy and hence a z = j z /B z , 
exactly for a force -free field 
(lAbramenko et al.ll 19961: iBao & Zhanglll998l) . 
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2. Compute a linear force-free extrapolation from B z and 
choose the overall value, 0^^, which best repro- 
duces the observed B x , B v distribution over the re- 
gion dPevtsov et al.ll 1991 lLongcope etaill 19981: IZhangl 
l2006h . 

The st udies by lHagino & Sakurail (120041) and iBurnette et al.1 
(2004) show that both techniques are generally consistent. 
The key result of these observations is a robust hemispheric 
rule whereby the average a value is negative in the north- 
ern hemisphere and positive in the southern hemisphere, al- 
though there is significant scatter including a mixture of signs 
of a within single active re gions. This hemisph eric pattern 
in a has also been found by lPevtsov et al.l (12001 1) who recon- 
structed the radial and toroidal components of the global mag- 
netic field under simplifying assumptions. 

A trans-equatorial sign change in helicity is supported 
by numerous proxy obs ervations such as Ha images of 
active region structure (iHald 19271). in situ heliospheric 
measurements (|Smith & Bieberl 1 1993b . differential rotation 
(Berger & Ruzmaikin 20 0^), and filament/p rominence mag- 
netic fields dRustl 11967b iMartin et"aIl fT994). Using newly- 
developed simulations of the global coronal evolution, we 
have recently been able to reproduce the filament hemispheric 
pattern including exceptions (with 96% agreement), in a 
comp arison with 109 observed filaments (lYeates et al.ll2007l 
2008). In this letter we describe the distribution of current 
helicity in a 30-month simulation, which we hope to compare 
with new magnetic observations from the SDO (NASA Solar 
Dynamics Observatory) mission. 

2. CORONAL MODEL 

Our simulations of the 3D coronal field evolution 
(I Yeates et al.l l2008) us e the coupled flux transport and mag- 
netofrictional model of|van Ballegooiien et al.1 (120001) . in a do- 
main extending from 0° to 360° in longitude, -80° to 80° in 
latitude, and R & to 2.5Rq in radius. The coronal magnetic 
field B = V x A evolves via the non-ideal induction equation 
dA 

— = vxB-r/ e j, (2) 
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FIG. 1. — Simulated magnetic field on days 100 (top) and 910 (bottom). 
Left column shows radial magnetic field on solar surface (white for positive, 
black for negative), and right column shows selected field lines of the 3D 
coronal magnetic field. 

in response to flux emergence and advection by large-scale 
motions on the photospheric boundary. Rather than solve the 
full MHD system we approxi mate the momentum equation by 
the magnetofrictional method (lYang et alH 986). setting 

1 i xB 

v=--^ + v out (r)r. (3) 

This artificial velocity ensures evolution through a sequence 
of near force-free states. The second term is a radial out- 
flow imposed only near to the upper boundary, where it sim- 
ulates the effect of t he solar wind in opening up field li nes in 
the radial direction (Mackay & van Ballegooiien 2006). The 
diffusivity r\ c consists of a uniform background term and 
an enhancement in regions of str ong current density j (see 
Mackay & van Ballegooiien 2006). 

The photospheri c boundary conditions are described in 
lYeates et al.l d2007l) : the surface flux transport model includes 
newly emerging magnetic bipoles based on active regions ob- 
served in synoptic normal-component magnetograms from 
NSO, Kitt Peak. The emerging bipoles take a simple math- 
ematical form, with properties chosen to match the location, 
size, tilt, and magnetic flux of the observed regions. They are 
inserted in 3D with a non-zero twist (magnetic helicity), cho- 
sen to match the observed sign of helicity in each hemisphere. 
The simulation illustrated in this letter models 30 months 
of continuous evolution during the rising phase of Cycle 23 
(from 1997 April 9 to 1999 October 10, rotations CR1921 
to CR1954). From an initial potential field extrapolation, the 
photospheric and coronal fields were evolved forward contin- 
uously for 914 days with 396 new bipoles inserted during this 
time. Two example snapshots of the simulated magnetic field 
are shown in Figure Q] 

3. SOURCES OF HELICITY IN SINGLE ACTIVE REGIONS 

To illustrate the sources of current helicity in our simula- 
tion within an individual active region, Figure|2]zooms in to a 
bipole in the northern hemisphere which emerged on day 136 
(as measured from the start of the simulation). 

There are three main sources of coronal currents and helic- 
ity in our model: 

1. The new bipoles emerge twisted. This twist is ini- 
tially concentrated low down in the centre of the bipole, 
as seen from the field lines in Figure |2ja) which are 




Lohgitude idcg) 

FIG. 2. — Structure of a single bipolar region, showing (a) magnetic field 
structure, (b) distribution of current helicity a on day 140, and (c) distribution 
of a on day 190. In (a) grey shading shows radial magnetic field strength on 
the solar surface (black negative, white positive), and coloured lines show 
selected coronal field lines. In (b) and (c) contours of a at height 14 Mm 
are shown in colour scale, and green contours show strength of radial surface 
magnetic field (solid for positive, dashed for negative). 

skewed as they cross the bipole's central polarity inver- 
sion line (PIL). The sigmoidal concentration of negative 
a at the centre of the bipole is clearly seen on day 140 
in Figure |2{b). 

2. When the bipoles emerge they displace older fields and 
produce currents at the interface b etween old and new 
flux systems (see lYeates et al.l2008l) . In Figure|2jb) this 
is visible at the NW edge of the new bipole where it ad- 
joins a pre-existing bipole, and a layer of positive a has 
developed. Note that this is opposite in sign to that from 
the twist of the new region, as seen in Figure[2fa). This 
corresponds to field lines that are oppositely skewed at 
this edge of the new bipole, as compared to those across 
the central PIL. This is just one example of how both 
signs of a may naturally be produced within a single 
active region, as found in observations. 

3. Over time, surface motions shear the coronal field gen- 
erating further currents. This is visible in Figure EJc), 
which shows the distribution of a for the same region 
on day 190, after 50 days evolution. There is a signif- 
icant build-up of negative a, particularly at the North 
and South ends of the bipole where helicity was initially 
low. This build-up is caused by differential rotation and 
convergence (due to supergranular diffusion). 

In addition to these sources of current helicity, it may also 
be locally reduced by diffusive cancellation and reconnection. 
Also, helicity is periodically removed through the top bound- 
ary of the domain when excessive build-up of twist leads to 
localised temporar y losses of equilibrium, and the eje ction of 
twisted flux ropes (Mackay & van Ballegooiien 2006). 
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FIG. 3. — Global distribution of a at height 14 Mm on days 10, 100, and 910. Left column shows contours of a (white for positive and black for negative, 
saturation level ±20 X 10~ 8 mT 1 ). Right column shows latitudinal profile, averaged over longitude in 2° latitude bins. Error bars show one standard deviation. 



4. GLOBAL DISTRIBUTION OF CURRENT HELICITY 

The global distribution of current helicity, a, is shown in 
Figure [3] at days 10, 100, and 910 of the simulation. From 
the initial potential field on day (with a = everywhere), a 
pattern of intermixed positive and negative a has developed 
by day 10, simply due to photospheric shearing — this is be- 
fore the first active region emergence. After about 100 days, 
a clear latitudinal trend in a emerges, although there is still 
significant local variation in both strength and sign. This pat- 
tern persists for the rest of the simulation, and up to medium 
heights in the 3D corona (nearer the top of the computational 
box high values of a become localized to closed field regions, 
with a k where the field is open). 

In Figures Oa), (e), and (f), it can be seen how the mean 
a at low latitudes (0° to about 50°) develops into the ob- 
served hemispheric trend, although with considerable scat- 
ter as observed on the real Sun. However, at high latitudes 



the sign of a is reversed. These polar reversals correspond 
to the East- West PILs at the polar crown boundaries, and 
move steadily poleward through the simulation as the polar 
crowns reduce in size towards polar field reversal (we are 
approaching solar maximum). This opposite sign of a is 
caused by differential rotation of the predominantly North- 
South field lines at this latitud e, and is a well-documented 
problem for theoretical mod els (Ivan Bal legooiien & Martens! 
1990; Rust & Kumadll994t) . At lower latitudes, as was il- 
lustrated by Figure Gfc), differential rotation of North-South 
PILs produces th e observed hemispheric sign of helicity 
ilZirkeretalJll997b . 

Figure [3] shows mean values of a at active latitudes of the 
order 10~ 7 m _1 . The actual maximum and minimum values 
recorded on day 910 of the simulation were 2.24 x 10~ 6 m _1 
and -1.84 x 10~ 6 m _1 . A key result of this study is that 
these values are much higher than those estimated from linear 
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force-free extrapolations. Such solutions suffer a constraint 
on the maximum a in order to obtain a decay with height 
dAulanier & Demoulinll 1998b . requiring that a < 2n/L x (the 
"first resonant value"), where L x is the horizontal length of the 
periodic box. The linear force -free model of an observed fila- 
ment by |Xulanieretaflj2(^ has a = 2.3 x 10~ 8 m _I , and for 
the solutions of iMackav et al.l (1 1999b this first resonant value 
was at a = 4.24 x 10 m _1 . By contrast, studies using nonlin- 
ear force-free extrapolations from vector magnetog rams us- 
ing the Grad-Rubin type method ( Ama riet al.ll 1 9971) find lo- 
cally higher values of a (e.g.. lBlevbel et al1l2002l) . They are 
also more realistic because they allow variable a within a sin- 
gle r egion, as in our simu lations. For a particular active re- 
gion, |R6gniereral] (2002) found maximum values of the or- 
der 10~ 6 m _1 , consistent with the results of our simulations. 

5. DISCUSSION 

In this letter we have shown how our 3D simulations of the 
global coronal magnetic field evolution are able to model the 
development and transport of current helicity, a, over many 
solar rotations. We find a clear latitudinal pattern of a that 
persists throughout the simulation, although locally within 
single bipoles there is significant scatter and intermixing of 
both signs of a, in agreement with observations. Local values 
may be much higher than those predicted by linear force-free 
extrapolations. 

With existing measurements of a limited to vector magne- 
tograms of individual active regions, robust observations of 
the latitudinal distribution of a await full-disk vector magne- 
tograms. These will shortly be available from the NASA Solar 
Dynamics Observatory (SDO) satellite. In particular the HMI 
(Helioseismic and Magnetic Imager) instrument will provide 
synoptic full-disk vector magnetograms at 1" resolution and 
approximately 90s cadence. This will offer an exciting oppor- 
tunity to test and refine our theoretical model for the coronal 
magnetic field. In particular, consistent measurements over a 



large portion of the solar cycle will allow us to consider how 
the helicity distribution varies over both space and time. 

Whether there is a systematic variation in the latitudinal 
tren d of helicity over the solar cycle remains a n unresolved is- 
sue (ISokoloff et al.ll2006tlPevtsov et al.ll2008|), and has impli- 
cation s for the sub-surface or i gin of helicity (IChoudhuri et al.l 
120041) . Indeed lKleeorin et all (12003) showed that observations 
of a in active regions provide importan t constraints o n theo- 
ries of the solar dynamo itself (see also Sokoloff 2007). Ejec- 
tion of helical fields from the corona, as included in our simu- 
lations, is also t hought to play an important role in sustaining 
the solar cycle (Blackman & Brandenburg 2003). 

A particular feature of our results is the sign reversal of cur- 
rent helicity at the high-latitude polar crowns. This would ap- 
pear to be in conflict with observations of magnetic fields in 
polar crown fila ments, which show no such reversal in their 
chirality pattern dRustj|1967t iLerov et all 1 19831; iMartin et all 
1994). We hope to address this outstanding issue in longer 
simulations covering a greater portion of the solar cycle. It 
is not at present clear whether longer-term poleward transport 
of the correct sign of helicity will be enough to counteract 
the effect of differential rotation on the North-South oriented 
field lines at these latitudes. Observations of vector magnetic 
field s in the polar reg ions, such as those being made by Hin- 
ode dLites et al.ll2008l) and soon the SDO mission, should help 
to constrain our models. 
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